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,— I Abstract. We consider the wave equation (9| — Ag)u{t,x) — f{t,x), 

in W^, w|]K_xK" = 0, where the metric g = (5'jfc(a;))"^_i is known outside 
an open and bounded set M C M" with smooth boundary dM. We define 
^ a deterministic source f{t,x) called the pseudorandom noise as a sum of 

^ point sources, f{t,x) = YlJLiC'j^xj{x)6{t), where the points xj, j G 

form a dense set on dM. We show that when the weights aj are chosen 
appropriately, mIrx^m determines the scattering relation on dM, that is, it 
^ determines for all geodesies which pass through M the travel times together 

^ with the entering and exit points and directions. The wave u{t,x) contains 

the singularities produced by all point sources, but when aj — X~^^ for some 
^ A > 1; we can trace back the point source that produced a given singularity in 

^ the data. This gives us the distance in (M."-,g) between a source point Xj and 

^ an arbitrary point y G dM. In particular, if {M ,g) is a simple Riemannian 

2 manifold and g is conformally Euclidian in M, these distances are known 

^ to determine the metric g in M. In the case when {M,g) is non-simple we 

• '-j present a more detailed analysis of the wave fronts yielding the scattering 

*%> relation on dM. 
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1 Introduction 

In this paper we consider an inverse problem for the wave equation 

{dl - Ag)u{t, x) = f(t, x) in (0, oo) x R", 
u\t=o^ dtu\t=o = 0, 

where Ag is the Laplacc-Bcltrami operator corresponding to a Riemannian 
metric g{x) = [9jk{x)]l^^^, that is 

j,k=i ^ ^ 

where \g\ = det{gjk) and = g{x)~^ is the inverse matrix of 

We assume that gjk G C°^(]R") and that there are Ci, C2 > such that 

n 

ciier < E 9A^)^'c' < C2ier, e M". (1) 

Moreover, we assume that the metric g is known outside an open and bounded 
set M C having a C°° smooth boundary dM. 

Denote by dM{x,y) = dM,g{x,y), x,y G M, the distance function of 
Riemannian manifold {M,g), where g is considered as its restriction to M. 
Let T > diam(M), where diam(Af) = max{(ijv/(a;, ?/); ,t, ?/ G M}. 

We choose the origin of the time axis so that the source / is active at 
time t — 0. To ensure compatibility with the the initial conditions we let 
To < and define the measurement map L — Lg, 

L:Cr(To,T)®Cr(M")^C°°((To,T)xaM), Lf = uyo,T)xaM: (2) 

where u is the solution of the wave equation 

{d^ - Ag)u{t, x) = fit, x) in (To, T) x M", (3) 

u\t=To = dtU\t=To = 0. 

Above, C^{To,T) denotes the space of smooth functions having compact 
support in {Tq,T). Its dual space, the space of generalized functions or 
distributions, is denoted by V^Tq, T). Moreover, for functions G C^(To, T) 
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and Tp G C^(]R") we denote their pointwise product by (0 ® ip){t,x) = 
(j){t)^/j{x). 

We remark that the assumption Q together with the finite speed of 
propagation for the wave equation imply that the measurement Lf does not 
depend on gjk{x), for |a;| > R, when R is sufficiently large. Thus we may 
assume without loss of generality that all the partial derivatives d^Qjk are 
bounded on M". 

Let Xj G dM, j = 1, 2, . . . , be a dense sequence of points in dM, and let 
us consider point sources 

f^^{t,x):=S{t)S,^{x), J = 1,2,.... 

In order to study the measurements Lf^j, we will use Sobolev spaces, see 
[56], 

H;{R') := {/ G S'iR'); Wf^.^^.^ := ||(1 - ^Y^' f\\,,^^,^ < +oo}, 
H;{U) :={feH;{R'y,snppfcU}, 
H;{U) := {/ G V'{U); f = h\u for some h G H;{R'')}, 

where f/ C M"' is open and s G M. When p = 2 we omit the subscript 
p in our notation, that is, we denote H^{U) = H2{U) etc. Moreover, we 
use projective topology on the tensor product X ®Y oi two Banach X and 
Y , that is, ||^||x(giy := i^^f X]j II II ^11%' II "^5 where infimum is taken over all 
representations z = '^jXj (8> yj. We also use projective topology on tensor 
products of locally convex spaces, see e.g. [23 Def. 43.2]. The measurement 
Lfxj can be defined in the sense of the following lemma. 

Lemma 1. Let p G (1, and let m ^ N satisfy m > Then the 

measurement operator L, defined in has a unique continuous extension 

L : H-\To,T) ® Hp\W) V\(Tq,T) x dM). 

We will prove Lemma [T] and other results presented in introduction in 
Sections IHMl 

In this paper we study a single measurement LHq that simultaneously 
combines all the measurements Lf^. by adding those together with appro- 
priate weights. When the measurements Lf^. are summed together, to the 
authors knowledge, there are no algorithms which could filter the value of 
a particular measurement from the sum. We will ask, however, can we find 
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the essential features given by these measurements, hke the travel times be- 
tween points on dM, so that the metric could be determined under certain 
geometric conditions. 

In Section[2| Definition[T| we construct a specific function ho{t,x), called 
pseudorandom noise, so that LHq determines the scattering relation T,M,g for 
the manifold (M, g). The scattering relation has been efficiently used to solve 
several geometric inverse problems [131 HHl EOl EI] . 

To define the scattering relation, let TM denote the tangent space of M 
and let 7 denote the tangent vector of a smooth curve 7 : [a, b] — )■ M. Let 
SM = {{x,^) G TM; W^Wg = 1} be the unit sphere bundle on M and define 

d±SM = e SM; X G dM, T{jy,Og > 0}, 

where u is the exterior normal vector of dM. Moreover, let t{x,^) be the 
infimum of the set {t G (0, 00]; 'yx,^(t) G dM}, where 'jx,^ denotes the geodesic 
with initial data (x,^) G TM. We define the infimum of empty set to be 
+00. 

The scattering relation is the map S = ^M,g, 
S : ^ d+SM X R, = {{x, G d^SM; t{x, < 00} 

defined by = (7x,5(t(x, 0), 7x,5(^(a;, 0), ^(^, 0)- 

Our main result is the following. 

Theorem 1. Let M C M", n > 2 be an open and bounded set having a 
C°° smooth boundary. Then there is a generalized function /io(t,x) such 
that it is supported on {0} x dM and has the following properties: As- 
sume that gjk,g'jk G C°^(]R") are two Riemannian metric tensors satisfy- 
ing Moreover, assume that gjkix) = g'jki^) for x G M" \ M. Let 
T > max(diam(M, (yf), diam(M, (7')), Tq < 0, and assume that 

Lgho = Lg^ho on (Tq, T) x dM. 

Then the scattering relations ^M,g one? T,M',g' of Riemannian manifolds (M, g) 
and (M, g') are the same. In particular, if (M, g) and (M, g') are simple, the 
restrictions of the distance functions on the boundary satisfy dM,g{x,y) = 
dM,g'{x, y) for x,y e dM. 

Recall that a compact Riemannian manifold (M, g) with boundary is 
simple if it is simply connected, any geodesic has no conjugate points and 
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dM is strictly convex with respect to the metric g. Any two points of a 
simple manifold can be joined by a unique geodesic. 

The key idea of proof of Theorem [l] is to use source /io(^, x) = Xljli (^jfxy 
The point source cij^f^.^ produces a smgularity, which is observed at a point 
y G \ M at time to = d{xjg,y) with a magnitude ajg/3{xj^,y), where 
P is an unknown nonvanishing smooth function. Appropriate choice of the 
weights ttj allows us find the index jo by looking nearby singularities. Indeed, 
when Xjf, — )■ xj^ and jk — oo, we see that the asymptotic behavior of the 
magnitude ajf,f3{xji^,y) as — )■ oo will be that of the weights Oj^.. Thus 
it is possible to factor out aj^ in the magnitude and determine aj^. This 
argument is presented in Section [7| and gives us the distances d{xj,y) in 
(M"", g) for arbitrary point ?/ G M" \ M and a source point Xj. 

Theorem [T] and boundary rigidity results for simple manifolds imply the 
following: 

Corollary 1. Let M C and let gj^^g'^^ E C°°(]R") be two Riemannian 
metric tensors satisfying assumptions of Theorem^ Assume that {M,g) and 
(M, g') are simple Riemannian manifolds. Then 

(i) If n = 2 and 

Lgho = Lg,ho on iTo,T) x dM (4) 

then there is a diffeomorphism $ : M — )■ M such that $|aA/ = Id and 
9 = '^*9'- 

(a) For n > 3 there is e = en,M > such that if \\gjk — 5jfc||c2(Af) < ^n, 
\\9'jk~ ^jk\\c^{M) < o,'>T'd ^) holds, then there is a diffeomorphism $ : M — )■ 
M such that ^\dM = Id and g = 

(Hi) If gjk{x) = a{x)6jk and g'^^ix) = a'{x)5jk, that is, the metric tensors 
are conformally Euclidian, and (01) holds, then gjk{x) = g'jki^) for x G M. 

Indeed, by Theorem [T| the case (i) follows from [l^, (ii) follows from [TU] , 
and (iii) from [Ml SQl HI] . 

If Uhlmann's conjecture [57], that the scattering relation determines the 
isometry type of non-trapping compact manifolds with non-empty boundary, 
can be proven, then Corollary [T] holds for more general class of manifolds. 

The problem of determining the metric g (possibly up to a diffeomor- 
phism) with given the measurement LHq with only one function ho{t,x) is 
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a formally determined inverse problem. Indeed, the formally computed "di- 
mension of the data", that is the dimension of (To,T) x 5M, is n and coin- 
cides with dimension of the set M on which the unknown functions gjk{x) 
are defined. 

The formally determined inverse problems have been studied in many 
cases. For instance, two dimensional Calderon's inverse problem [31 IH 1221 1131 
l53] is formally determined. The same is true for the related inverse problem 
for the Schrodinger equation in two dimensions [9] . The corresponding inverse 
problems in dimension n > 3, see [HI [281 EH 1121 El] and references in [16], 
are over-determined, that is, the dimension of the data is larger than the 
dimension of the unknown object. Similar classification holds for the elliptic 
inverse problems on Riemannian manifolds [171 [TBI IHH [35l |37| . Moreover, the 
inverse travel time problems, i.e. boundary rigidity problem, see [301 ES EH 
HHl HH HHl |52] , is formally determined in dimension n = 2 and overdetermined 
for n > 3. 

Inverse problems in time domain related to the Laplace-Beltrami operator 
Ag, namely the inverse boundary value problem for the wave, heat, and 
the dynamical Schrodinger equations with Dirichlet-to-Neumann as data, 
see [21 [71 [221 126], are overdetermined in dimensions n > 2. However, these 
problems are equivalent to the inverse boundary spectral problem, see [27] . 
and assuming that the eigenvalues are simple, Dirichlet-to-Neumann map at 
a generic Dirichlet boundary value determines the boundary spectral data 
[32l l33l Wl\ . Thus under generic conditions on the spectrum and on the 
boundary value (that is, under conditions that the these data belong in some 
open and dense set) it is possible to solve a formally determined inverse 
problem in time domain. 

We point out that in this paper we do not impose any generic conditions 
on the geometry and we give an explicit constuction of the boundary source. 
The boundary source considered in this paper is based on the idea of imitating 
a realization of white noise, and due to the many useful properties of white 
noise process, we hope that the constructed source may be useful in the study 
of other inverse problems requiring generic assumptions on the source. 

Another formally determined hyperbolic inverse problem, namely mea- 
suring Neumann data when the initial data {u\t=o., dtu\t=o) is non-zero and 
satisfies subharmonicity or positivity conditions, has been studied using Car- 
leman estimates [H [231 |29]. The present paper is closely related to these 
studies, but we emphasize that we assume that the initial data for u van- 
ishes. 
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2 Pseudorandom noise as a source 



In this section we define a special source ho{t,x) which we call the pseu- 
dorandom noise. The specific assumptions on the amplitudes are explained 
in Section [7j An important feature of the pseudorandom noise is that it is 
supported only on a single point in time. 

Definition 1. Let xj G dM, j = 1,2, ... , be a dense sequence of disjoint 
points in dM , and let aj G M, j = 1, 2, . . . , J^JLi Wj\ < oo be a sequence of 
disjoint numbers. 

We define the pseudorandom noise on {xj)'^i C dM with coefficients 
{aj)'^i C M as the following generalized function on M x M".- 

oo 

i=i 

where 6{t) and6xj{x) are Dirac delta distributions on M andW^, respectively. 

It is rather straightforward to show that Hq is well-defined. First, it is well 
known that 5{t) G H'^R) and6^^{x) G C(M")'. Moreover, H^,{W) C C{W) 
when 1 > n/p' due to [Ml Thm. 2.8.1]. According to pEl Thm. 2.6.1] the 
dual space satisfies {H^,{R'')y = H'^W) with 1/p' = l-l/p and hence 
C{W)' C Hj;^{W) for 1 < p < Since \aj\ < oo we have 

oo 

i=i 

This yields that for any p G (1, and e > the pseudorandom noise /iq 
satisfies 

hoeH-\-t,e)®H;\M). (5) 

The spatial structure of the pseudorandom noise can be motivated by the 
structure the white noise. In the 1-dimensional radar imaging models, white 
noise signals are considered to be optimal sources when imaging a stationary 
scatterer [19] . This is due to the fact that different translations of the white 
noise signal are uncorrelated. In a similar fashion we have the following 
property for the pseudorandom noise Hq: for each and each sequence 
(xj^,)^^ converging to and satisfying Xj^, ^ Xj^, for all k G 2,^, it holds 
that aj^ — i- 0. This property will be crucial in what follows. 



7 



Moreover, a natural strategy to choose the points Xj is by random sam- 
phng. The term pseudorandom refers to the fact that the algorithmic gener- 
ators of random numbers use, in fact, a deterministic function to produce a 
sequences of numbers through so mixing process, that the user of the algo- 
rithm can consider the numbers to be analogous to independent samples of 
a random variable. In this manner, the pseudorandom noise can be seen as 
an imitation of a realization of a noise process. 

Another source of inspiration for us was a rather new measurement paradigm 
called compressed sensing fl2[ [H] , where one aims for a sparse reconstruc- 
tions of a linear problem using a small number of noisy measurements. We 
point out that by using the pseudorandom noise one can compress the mea- 
surements Lfx^ with point sources fx^ into a single measurement Lho. 

3 Measurement map 

In this section we prove that the measurement is well-defined when we have 
pseudorandom noise as source. 

Next, we consider the operator W : f ^ u mapping / to the solution of 
the equation ([s]). We call such operator the solution operator for the equation 
g. First, we note that by [211 Thm. 23.2.2], the operator W : f ^ u 
mapping / extends in a unique way to a continuous linear operator 

W : L\{To,Ty,H'{W)) ^ C{[To,T];H'+\R'')), s G M. (6) 

Moreover, if / G C°°([To,T] x M") and supp(/) CC (To,T] x M", that is, 
supp(/) is a compact subset of (Tq^T] x M", then Wf G C°°{[To,T] x M"). 

We will compose the operator W with the one-sided inverse X of the 
derivative dt, which is given by 



One sees easily that this operator has a unique continuous linear extension 



Next we prove that the measurement map L has unique continuous ex- 
tension 




(To,T). 



H-\T^,T) ® H'\W) V\{To,T) x dM). 



(7) 
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Proof of Lemma^ For sufficiently large z G M+, operator z — Ag is an iso- 
morphism between spaces if*"'"^(M") and if*(M"') as well as between spaces 
H^+^{W) and H^{W') for all integers s by gH]. 

By the definition of L, we have that L = Tr o W, where Tr is the trace 
operator 

Tr{u) = n|(To,T)xaA/, u e C°°((To,T) x M"). 

Let / G C^{{To,T) x M"). Then the solution u = Wf oi the wave 
equation [d^ — Ag)u = f can be written in the form 

m—l 

Wf = {z- d^riz - Ag)--Wf + dly{z - Ag)-'-^f. (8) 

i=o 

Now / = (9iX/, where Xf is C°°-smooth and satisfies supp(X/) CC {Tq,T] x 
M". By g, VTX/ is C°°-smooth and dtWXf = WdtXf = Wf. Hence 

m—l 

Lf = dt{z - d^rTiiz - Ag)-^WXf + - d^YM^ - ^.jY'-'f- (9) 

i=o 

Let us next consider terms appearing in ([9]). First we consider extension 
of the operator 

N N 



Y^^k^^k^ _ d^YTiiz - Ag)-'-^{(Pk ® ^fc) (10) 

A;=l 
AT 

- dfy<P,) ® {Tt{z - Ag)-'-^^,)), J = 0, . . . , m - 1, 



fe=l A;=l 
AT 



k=l 



mapping C~(To, T) ® C~(M") to C°^((To, T) x ^M). By |56l Thm 4.7.1] the 
maps 

are continuous, where Bp^p^^^{dM) is the Besov space on DM. Thus the 
operator (10) has a continuous extension in spaces ([T]). 
Next, consider the extensions of the operator 

AT AT 

J] 0fc ® V-fc ^ 5] dt{z - dtri:i{z - Ag)-'^w{{X(t>k) ® i^k) (11) 

k=l k=l 
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mapping C~ (To, T) ® C~ (M") to (To, T)x9M). As > -1-nwe 
have by [56!, Thm. 2.8.1] a continuous embedding H-^{W) ^ H-^-''/^{R''). 

Moreover, the operator X : H~^{Tq,T) — > L'^{Tq,T) and the embedding 
L2(To,T) ® i7-i-"/2(^M") L2((To,T);if-i-"/2(^n)) g^^g continuous. Thus, 

by (@, 

VTX: #~^(To,T) ® i7-^(M") ^ C([To, T]; //-"/^(M")) 
is continuous. 

Thus, as (1 - Ag)-"" : C([To, T]; ff-"/2(M")) ^ C([To, T]; /f-n/2+2m^^n^^ 
is continuous and — n/2 + 2m > 1/2, we see that the operator 

Tr(l - Agy'^WI : H-\To,T) ® Hj;\W) ^ C([To, T]; ^^(aM)) 

is continuous. 

Combining the above resuhs, we see that the operator (|9| has a continuous 
extension to spaces ([T]). As the spaces C^(To,T) and C^{M."') are dense in 
H~^(Tq,T) and H~^{R"-), respectively, we see that the continuous extension 
of L is unique. □ 



4 Inner product of a solution and a source 

Lemma 2. Let f G C;;«((To,T) x M), to G (To,T) and let w G C°°([To,to] X 
M") satisfy 

[d'^t - Ag)w = 0, zn (To, to) x M". 



Then 

"to 
To 



fit, x)w{t, x)dtdV{x) 

{{dtWf){to,x)w{to,x) - {Wf){to,xMw){to,x))dV{x) 



where dV{x) = \g\^^'^dx is the Riemannian volume measure of (W'^,g) and 
W : f i-^ u is the solution operator of wave equation 
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Proof. By finite speed of propagation of waves, see e.g. [211 pp. 150-156], 
supp{Wf{t)) is compact in M". The claim follows by integration by parts 

{{dtu){to, x)w(to, x) - u{to, x){dtw){to, x))dV{x) 

{{dtu)(To,x)w(To,x) - u(To,x){dtw)(To,x))dV{x) 

{{d1 - Ag)u(t, x) w{t, x) - u{t, x) {d^ - Ag)w{t, x))dtdV{x) 
fit, x)w(t, x)dtdV{x). 



f 

(7b,io)x 



□ 



Next, we will prove a generalization of the previous lemma for non-smooth 
sources /. Denote by B{0,R) = {x E M"; |x| < R} the Euchdean ball. The 
finite speed of propagation for wave equation, yields that there is i? > such 
that all / e C~((To,T) x M) satisfy supp(iy/) CC (To,T] x B{0,R). We 
define 

n:= B{0,R)\M. (12) 

Below, we use the fact (see [151 Thm. 7.2.3/6, Thm. 5.6.3/6]) that the 
operator Wq : h v mapping h to the solution of the equation 

(d'^ - Ag)v{t, x)=0 in (Tq, T) x fi, (13) 

v\{To,T)xdn = h, 

v\t=To = 0, dtV\t=To = 0. 

is continuous Wn : C^{{To,T) x dn) C°^{[To,T] x H). 
We let to £ (^0) T) and denote 

S := {to} X n (14) 

We denote the trace on S by Tr^, that is, we define {Tt-eu){x) := u(to,x). 
Let u = v{z) denote the exterior unit normal vector of dM at z. 

Moreover, let U be an open subset (or a submanifold) of M" and let us 
denote by dV (or dS) the Riemannian volume measure of {U,g). We embed 
the test functions into the spaces of distribution by using the inner product 
of the space L'^iJJ; dV), that is, we identify u G C^{U) with the distribution 

i> ^ u{x)iIj{x) dV{x). (15) 
Ju 
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We will denote the distribution pairing of u G V'{U) and ip G C^{U) by 
(u, ip)v'{u) and use analogous notations for other distribution pairings. 



Lemma 3. Let tg ^ {Tq,T) and define S by (14)- Then operators Tty^Wq 



and TiY;dtWQ have unique continuous extensions £'({To,to) x Oil) V(^l). 



Proof. Let v satisfy (13). Consider a function w G C°°([To, to] x ^l) such that 



{dt - Ag)w = in (To,to) x ^ and M;|(To,to)xan = 0. Then 



= / {{dt - Ag)v)w - v{{dt - Ag)w)dV{x)dt 

t=to p 

+ / {{duv)w — v{duw))dS{x)dt 

t=To JdQx{To,to) 



I 

Jn 



{{dtv)w — v{dtw))dV{x 
{{dtv)w — v{dtw))dV{x) 



h{dyw)dS{x)dt, 

t=to JdQx{To,to) 



in 

where d^, is the normal derivative on dQ, 

Denote by Wi : fi ^-^ w the solution operator of the equation 

{dt - Ag)w{t, x) = in (To, to) x fi, 

1^\{To,to)xdn = 0, 

w\t=to = fi, dtw\t=to = 0. 

Operator Wi : C^{Q) — )■ C°°{[To, to] xQ),is continuous, as can be seen using 
[T5| Thm. 7.2.3/6, Thm. 5.6.3/6]. Hence the operator 

d,W, : C^{Q) ^ C°°([To,to] X dQ), f ^ dyWJU 
is continuous. Moreover, 

L'^{n;dV) — 

{h,d,WJ\) 

L^{(To,to)xdn;dt(^dS)y 

where du is the normal derivative on dQ. We define the extension of TTj^dtW^ 
by identifying it with the transpose {d^WiY : £'{{To,to) x dQ)) V'{VL) of 
the operator d^Wi : C^{Q) C°°([To,to] x dQ). 

Similarly, we define the extension of Ti^Wn by the transpose {dyW2y : 
£'((To,to) X dQ)) V'{Q) of dyW2 : C^{Q) -> C^{[ToM x dQ), where 

: /2 ^ w is the solution operator of the equation 

{dt~Ag)w{t,x) = Q in(To,to)xfi, 

w\{TQ,to)>^dQ. = 0, 

w\t=to = 0, dtW\t=to = -f2- 

□ 
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Denote by d^{x,y), x,y & Q, the distance function of Riemannian mani- 
fold (0,5^1^). Next we generalize the result of Lemma [2] for a larger class of 
functions. 

Lemma 4. Let to G (0,T) and e > satisfy [— e, e] C (To, to)- Define E by 
Q. Let f e H-\-e,e) O H~\M) and w E C°°([To,to] x M'^) satisfy 

{dl - Ag)w = 0, zn (To, to) x M". 

Suppose that w{to),dtw{to) G C^(fi), anc? /ei x ^ C'^(^05^o) satisfy x = 1 
m a neighborhood of [—e,to — r], where 

r := dn[supp{w{to)) U supp(9tw;(to)), f^^)- 

T/ien 

(/, w)£-'(R"x(To,io)) (16) 

= iTT^dtWnxLf,w)T)'(Q) - (TrsV^nx^/, 5tW^)©'(n), 

where we have defined Lf = on dB{0, R). Here we regard Q as Riemannian 
manifold {Q, g\n)- 

Proof. We suppose first that / G C*^((— e, e) x M). Recall that W is solution 
operator of wave equation Then Wf{-,t) = if t < — e, and 

Lf = TianWf = x^ranWf, in (Tq, to - r) x dQ, 

where Tian is the trace on (To,T) x dQ. As Q (1 M = 0, we have that 
(92 - Ag)Wf = in (To, T) x Q. By uniqueness of the solution of ^ 

Wnx'TranWf = Wf, in (To, to - r) x Q. 

By finite speed of propagation 

Ti^dlWnx'^idnWf = lij^dlWf, j = 0, 1, 

on {to} X supp(w(to)) U supp(9(w(to)). By Lemma[2| (16) holds. 
Then the claim follows as the embeddings 

are dense and operators {Trj:diWn)xL : H~\-e, e)^H-^{M) V'{{To,to)x 
dQ), j = 0,1, are continuous. □ 
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5 Gaussian beams 



We consider solutions of wave equation which are known as Gaussian beams 
[SI El US] • These solutions have been constructed to analyze the propagation 
of singularities for the wave equation in the presence of caustics. Here we 
use Gaussian beams as an auxiliary technical tool to analyze singularities in 
the measurements. 

Definition 2. Let e > 0, G N and let 'j be a unit speed geodesic on 
(M", g). A formal Gaussian beam of order N propagating along geodesic 7 is 
a function of form 

N 

U^{t,x) = e-"/^exp{-(ze)-i^(t,x)}^M„(t,x)(ie)™, t e M,x G M'^ 

m=0 

satisfying the following properties: The phase function 9 and the amplitude 
functions Um, rn = 0,1, . . . , N , are complex valued smooth functions. The 
phase function 9 satisfies the conditions 

9{t, 7(t)) = 0, Im 9{t, x) > Co{t)d{x, ^{t)f 

where Co(t) is a continuous strictly positive function. The amplitude function 
uq satisfies uo(t,'y{t)) 7^ 0. Finally, for any compact set K CC M x M" there 
is a constant C > such that the inequality 

\{d^ -Ag)U^{t,x)\ <Ce^-'^/^ 

is satisfied uniformly for {t, x) G K. 

The construction of a formal Gaussian beam U^{t,x) is considered in 
detail e.g. in Sect. 2.4]. Next, we recall the construction and pay attention 
to the properties of Gaussian beams which we need later. 

Let us write the geodesic 7 in the usual coordinates of as 7(t) = 
(7^(t), . . . , 7"(t)). We construct the phase function 9{t, x) at each time t G M 
in terms of a finite Taylor expansion in the x variable centered at 7(t), 

9{t,x)= J2 ^(^-7wr, 

\a\<N 

where 9a are smooth functions and G N. 
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Let Cj = {Sij, . . . , Snj) be multi-indexes with the value 1 at the jth place. 
For clarity, we use the notation pj (t) = 6^^ (t) for the first order coefficients 
and the notation Hjk{t) = 6a{t), a = ej + e^, for the second order coefficients 
in the expansion of 6. 

The construction of a formal Gaussian beam consists of the following 
steps. 

1. We define Ooit) = and pj{t) = ELi ^ifc(7(^))7''(^)' ^hat is, the first 
order coefficients Pj{t) are the covariant representation of the velocity 
vector 7. 

2. The symmetric matrix H{t) = [-f^jfc(^)]"fc=i of the second order coef- 
ficients are obtained by solving a Riccati equation, or an equivalent 
system of ordinary differential equation. We write H{t) = Z{t)Y{t)~^ , 
where the pair of complex nx n matrices {Z{t), Y{t)) is the solution of 
the system of ordinary differential equations, 

jY{t) = BityVit) + C{t)Z{t), Y\,^, = Y\ 

j^Z{t) = -D{t)Y{t) - B{t)Z{t\ Z\,^, = Z'. 

Here we choose the initial values to be Z^ = il and Y^ = I, where / 
is the identity matrix and i is the imaginary unit. The matrices B{t), 
C{t), and D{t) in M"^" have components given by the second deriva- 
tives of the Hamiltonian h{x,p) = k=i9jki^)P'p''y^'^ evaluated in 
the point {x,p) = {■y{t),p{t)) 

^ dx^dpj ' dpjdpi ' dxWx'- 

The fact that the complex matrix Y{t) is invertible for all t G M is 
crucial for the construction, and is discussed in detail in [261 Section 
2.4]. 

3. The coefficients 6a{t) of order |a| = m > 3 are solved inductively, 
with respect to m. The coefficients 6a{t) are constructed using the 
coefficients Oa{t) defined so that 

\a\=m \a\=m 
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for all y = Y-'^{t){x - j{t)), y e C". We obtain the coefficients Bait) 
by solving successive linear systems of ordinary differential equations 

j^ea{t) = K^{t), ea{o) = o 

where Ka{t), depend on 9p{t) with \(3\ < m — 1, the matrix H{t), vector 
p{t), and the metric gjk and its derivatives at "y{t). 

4. When the phase function 9{t, x) is constructed, the amplitude functions 
Un{t,x) are solved using the transport equations, or equivalently, the 
following ordinary differential equations. Let 

\a\<N 

where the coefficients Um,ait) are obtained by solving the successive 
equations 

where r{t) and J^m^aif) depend on Um',p with \[3\ < \a\ +2 and m' < 
m — 1, the function 9{t, x), the metric Qj}^ and their derivatives at (i, x), 
X = 7(i). 

By the above construction, we have the following remark. 

Remark 1. The phase function 9{t,x) and the amplitude functions Um{t, x) 
at time t — have the form 

n 

9{0, x)^J2 93k{y)r]^{x^ - y') +i\x- y|^ 
j,k=i 

where {y,r]) = (7(0), 7(0)) is the initial data of the geodesic 7, uq{0,x) ~ 1, 
and Um{0,x) = for m > 0. Hence U^{0,x) is dependent on the metric 
Qjk only via gjk(y). Moreover, dtU^ {0,x), although of more complex form, is 
dependent on the metric gjk only via d°'gjk{y) for a certain finite collection 
of multi-indices a G N". 
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If the coefficients of an ordinary differential equation depend smoothly on 
some parameter so does the solution P, and thus we see using an induction 
that the phase function 9 and the amplitude functions Um depend smoothly 
on the initial data (y,//) = (7(0), 7(0)) of the geodesic 7. In particular, the 
amplitude function Mo(t, x; y, r]) satisfies 

uo G X M" X SW). (17) 

To this far we have considered a formal Gaussian beam. By using con- 
tinuous dependency of the solution of the wave equation on the source term, 
on obtains the following results, see e.g. [2S]: 

Let 7 be a unit speed geodesic, G N, e > and let be a formal 
Gaussian beam of order propagating along geodesic 7. Let x ^ (^^(M") 
be a function which is identically one in a neighborhood of 7(0) and let to > 



and let R be the radius in the equation (12). Then for j G N and a G N" 
satisfying j + |a|<A^ — r;,/4 there is C > such that the solution of the 
wave equation 

[d^ - Ag)w,{t, x) = 0, (t, x) G (To, to) X M", (18) 
Weito,x) = xix)U^{0,x), dtW,{to,x) = -xix)dtU^{0,x). 

satisfies 

sup \d{d^{w,{to -t,x)- f/f (t,a;))| < Ce^-O+I"!)— /4. (19) 
xeB{o,R),te(To,to) 

We call a Gaussian beam of order propagating along geodesic 7 back- 
wards on time interval (To, to). 



6 Determination of the travel times 

Lemma 5. Let be a Gaussian beam of order > 1 + n/4 propagating 
along geodesic 7 backwards on time interval (To, to), that is, let be the 
solution of (18). Let Hq be the pseudorandom noise 

00 

ho{t,x) = Y,(^jHt)^x,{x). (20) 
i=i 

If 7(^0) 7^ Xj for all j = 1,2, . . . then 

lime''/^(/io, W£)£-'{i;"x(ToA))) = 0- 



17 



Moreover, if •yito) = Xj then 

linie"/^(/io,We)£-'(Rnx(To,to)) = ajUoito,Xj)\g\^^'^ixj), 

where UQ{t, x) is the first amplitude function of a formal Gaussian beam prop- 
agating along geodesic 7. 

We remind the reader that the test functions are embedded in S'{W^ x 
(To,T)) using (pi). 



Proof. By equation ( 19 ) we have that 



00 

= ^OjUoito, Xj) exp {-{iey^e{to, Xj)}\g\^/'^{xj) +0(e). 
i=i 

As lm9(to,Xj) > Co(to)d{xj,'y(to)) we have that 

\exp{-{iey^9{to,Xj)}\ = 0(e), if 7(^0) 7^ x^. 

Suppose that 7(^0) = ^j- Then exp { — {ie)~^9(tQ,Xj)} = 1 and there is a 
constant C > depending on 7 and to such that 

\e"'^'^{ho,We)s'iR"x(To,to)) - ajUo{to,Xj)\g\^^^{xj)\ 

<CY,Mexp{-{te)-^e{to,Xk)}\+C J2 \ak\\exp{-{te)-'e{to,Xk)}\ 

k=l k=j+l 



+ C^|ai| +0(e). 



We may first choose large / G N and then small e > so that the above 
three sums are arbitrary small. The case, 7(^0) 7^ for all j = 1,2,..., is 
similar. □ 

Next we define an auxiliary function S{yQ, tjq, tg) which is non-zero if and 
only if there is j G Z+ such that 'yyo,r]o(to) = ^j- 
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Definition 3. Let {yo,r]Q) E TM" be such that yo E and \\r]o\\g = 1. 
We denote by 'y{t]yo,rio) = 'jyo^noit) the geodesic on (W^,g) with 7(0) = yo, 
7(0) = rjo. Moreover, let be a Gaussian beam of order N > 1 + n/4 
propagating along 'y{t;y,ri) backwards on time interval (To, to)- We define 



Lemma 6. Let {yo,r]) E SQ and to E (0,T). Then Lho for pseudorandom 
noise ho and {Q,g\Q), given as a Riemannian manifold determine S{yo, rjo, to). 

Proof. Let be a Gaussian beam of order N > l+n/4 propagating along tlie 
geodesic 7(-;i/0)^o) backwards on time interval (To,to)- We may choose the 
cut-off function x in the equation (fl8| so that We{to), dtWe{to) E C^{Q). As 
g\fi is known, we have by Remark 111 that the initial data We{to), dtWe{to) are 
known. Moreover, operators Ti-sd^Wn, j = 0, 1, E := {to} x Q, are known. 
After choosing a suitable cut-off function x in Lemma |4] we have that the 
measurement Lho determines the distributional pairing (/io) 'W^e)£■'(M"x(To,^))• 
Hence S'(?/o, ?7o, to) is determined. □ 

The implicit function theorem yields the following remark. Note that 
to € M in the remark is not necessarily the first intersection time. 

Remark 2. Let {yo,f]o) ^ *S']R" and to G M satisfy 



Then there are neighborhoods / C M and U C SW^ of to and {yo,f]o) o-nd a 
smooth map i : U ^ I such that for t E I and {y, t]) E U 



We remind the reader that t{x,^), {x,C,) E TM", is defined as the first 
intersection time with dM, that is 



In the following, we use the Sasaki metric on the tangent bundle TM. 



S{yo,Vo,to) ■■= lime^/^l/io, w^£)£-'(R"x(To,io)) 



h{to;yo,Vo),i{to;yo,Vo)) e d±SM. 




r{yo,Vo) ■= inf{t E (0, 00]; 7(t; yo, ^0) e dM}. 
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Lemma 7. The first intersection times t : SQ — )■ (0, oo] and r : d^SM — )■ 
(0, oo] are lower semi- continuous. 

Proof. Let us consider r on SQ. Let a sequence {{yj,Tij))'^-^ C SQ converge 
to {yo, rjo) G SQ as j — )■ oo. We denote 7j(t) := 7(t; yj, rjj) and Tj := r(?/j, r/j). 

We will show next that lim inf j_i.oo ^ (0,'ro). Let t G (0,ro). Then 
70 (t) ^ and 

d(7o(t),5M) >0. 

Let j G Z+. Suppose for a moment that < oo. Noting that •jj is unit 
speed and "ijijj) G (9M, we have 

|t-r,| > rf(7,(t),7,(r,)) > rf(7,W,5M). 

If Tj = OO, then |t — Tj | = oo > d{'^/j{t), DM). 

The convergence 7j(t) — 7o(i), as j — )■ oo, implies that for large j 

\t - Tj\ > rf(7o(t),(9M)/2 > 0. 

Hence liminfj_i.oo tj ^ t for all t G (0, tq). 

Clearly liminfj_j.oo Tj > 0, and there is J G Z+ such that 

Tj > d{yj, dM) > d{yo, dM)/2 > 0, j > J. 

Hence liminfj^oo and liminfj_i.oo tj > tq. 

Let us consider r on d^SM. Let a sequence {{yj,rij))'jLi C d^SM con- 
verge to (|/o,?7o) G d-SM as j ^ oo. We denote 7j(t) := 'y{t]yj,rij) and 
:= r(?/j-,r/^). 

Repeating the above argument, we see that liminfj_^oo Tj ^ (0, tq). Thus 
it is enough to show that liminfj_>.oo tj ^ 0. 

Remark [2] gives neighborhoods / C M and U C SW!^ of zero and {yo,i]Q) 
and a map i : U — )■ J of boundary intersection times. We denote V := 
t/ n d^SM. As 7(0; x,0 e dM for e 1/, we have £ = in V". In 

particular r := d{i{V),R\ I) > 0. For large j, (7j(0), 7j(0)) G V and thus 

7i(t)eM, tG(0,r). 

Hence tj > r > for large j, and liminfj^oo > tq. □ 
We easily see the following continuity result for r. 
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Figure 1: On left, trajectory of a Gaussian beam propagating along geodesic 
7(t) := 'y(t;yj,rij) backwards on time interval (To,tj). If S{yj,rij,tj) 7^ 0, 
then there is a point source at •yitj). On right, a sequnce {yj,rij) e Sfl 
converging to (x, ^) G d-SM and trajectories of the corresponding geodesies. 



Lemma 8. Let {{yj^rij))"^^ C SVL converge to (x,,^) G d^SM in the Sasaki 
metric. Then lim^^oo Vj) = 0- 

Theorem 2. Let (a;,^) G d-SM and denote by J{x,^) the set of sequences 
{{tj]yj,i]j))'^i C (0,00) X SQ for which 

\im{yj,T]j) = lim tj G (0,oo), S{yj,rij,tj) ^ 0. 

The function S : Sfl x (0, 00) — )■ C determines t : d-SM — > (0, 00] 6t/ ^/ie 

r(x,0 = inf{hm tj; {{tj;yj,Vj))T=i ^ ^i^^O for some {{yj,Vj))T=i ^ 3^}- 

Moreover, if t{x,^) < 00, then there is a sequence {,{tj]yj,Vj))^=i ^ 
J(x,^) satisfying 

j-s>oo 

Proof. Let (x,^) G d_SM and {(tj;yj,rij))'jL-^ G J(x, ^). Let us show, that 
T(a;,^) < lim j^^tj. By Lemma |8| := T(yj,r]j) — )■ as j — )■ 00. We define 
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As S{yj,rij,tj) 7^ 0, we have 

7(^i - Tj-; Vj, Q = litf, Vj, Vj) e dM. 

As limj_^oo tj > and \imj_^oo Tj = 0, we have tj — tj > for large j. Thus 
0) — ~ '0' f*^^ large j. Moreover, 

lim = (7(0;a;,e),7(0;x,O) = 

In particular, {yj,^j) E d-SM for large j. Hence Lemma [T] gives 
lim tj = lim {tj — Tj) > liminf T{yj, ^j) > r(x, ^). 

In particular, we have proved the claim in the case t(x,^) = 00. 

Let us assume that t{x,S,) < 00. It is enough to show, that there is a 
sequence {(tj;yj,rij))'jL^ G J{x,C,) satisfying r(x,^) = limj^ootj. We denote 

to:=r(a;,0, z := -f{to;x,^), C := -7(^0; a;, 0- 

We have 

{x, = (7(^0; z, C), -7(^0; z, 0). 
As (rc,^) G d^SM, Remark |2] gives neighborhoods / and U of to and (zX) 



and a map £:[/—)■/ of boundary intersection times. After choosing local 
coordinates around z we may define 

where {xk^ ) C t/ is a subsequence of the dense sequence of source points 
in (20) satisfying limj^aoXkj = z and (tj)jli C / satisfies 

tj > i{xk^,C), lim tj = i{z, C) = to- 



Clearly {{tj;yj,r]j))°°^^ G J(x,0 and 



lim t,- = to = r(x, ^). 

j->oo 



□ 
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7 Determination of the scattering relation 

In the next theorem we consider pseudorandom noise hQ{t, x) with coefficients 

= X-^\ 

with some A > 1 and make computations "modulo an error in A" , where 

A = {-X^ : J G N}. 

For this end, let 771^(3) be the real number r such that s = r + a where a E A 
and r has the smallest possible absolute value. In the case when both r and 
— r satisfy this condition, we choose the positive value. 

Lemma 9. Let {yo,f]o) ^ SQ, to ^ (0,7"), and suppose that S{yo,riQ,to) 7^ 0. 
Then there is a sequence {{yj.rjj))^^ C 5*^2 and {tj)f=i C (0,T) such that 

iVj^Vj) ^ iyo,Vo), tj^h, S{yj,rij,tj) ^ 0, as j 00, (21) 

Suppose, moreover, that the coefficients of the pseudorandom noise ho are 
aj = A~'^\ Then for any sequences {{yj,rij))'^i C TM" and (tj)jli C (0,T) 
satisfying (21) we have that 

lim m^(log;^ \S{yj,T]j,tj)\) = log;, |?^o(^o, 7(^0); l/o, ^o)|^r^^(7(^o))|, 



where 7(t) = 7(t;yo)''?o) and uq is defined as in (17). 
Proof. We will use notation 

7i(*) := 7{t;yj,Vj), Zj := -fj{tj), Sj := S{yj,T]j,tj), 
(3j := \uo{tj,zfyj,r]j)\g\^'^{zj)\. 

As 5*0 7^ 0, we have that Zq = Xj for some j = 1, 2, . . . . By continuity of 
the geodesic flow and and density of (xj)'^^ C dM, there exists a subsequence 
C (a;,)~i and sequences iiyj,Vj))T=i ^ ™" and C (0,T) 

such that 

Xkj zo, {yj, Tjj) {yo, 770), tj -> to, as j ^ 00 
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and Zj = Xk^ ^ zq. Then \Sj\ = |afcj/3j ^ 0. As ^ zq and Xk^ zq, we 
have that kj — ?■ oo and thus ak — ?■ 0. By (17) and continuity of the geodesic 
flow, it holds that (3j — )■ /3o > 0. Hence Sj — )■ 0. 

Next we use the assumption that aj = X^^\ Let {{yj,rij))'^^ C TM" 
and {ij)^=i C (0,T) satisfy (21). As 5*^ 7^ we have that \Sj\ = ak^Pj for 
some subsequence (0^^)°^! (1(0^)°^]^. As Sj — )■ we have that ak^ — )■ 0. 
Moreover, sequence (logg Pj)'jLi is bounded. This boundedness together with 
logA«fcj e A and log;^aA;^. -00 yield 

mA(log^ ak^ + log;, (3j) = \og^ (3j 

for large j G N. Hence, 

hm mA(log;, \Sj\) = hm log^ l3j = \og^l3o. 

□ 

Theorem 3. // the coefficients of the pseudorandom noise ho are aj = X^^\ 
then the functions S : SQ x (0, 00) — )■ C and r : d^SM — )■ (0, 00] determine 
D(S) and 

Proof. Clearly r on d-SM determines -D(S). Let {x,C,) € -D(S). By Theo- 
rem[2|we may choose {{tj;yj,r]j))'^i € J{x,^) such that limj_s.ooij = ^). 
As 5\yj,7]j,tj) 7^ 0, l{tj]yj,r]j) = Xk^ for some subsequence (xfcj^^i of the 
sequence of source points. By Lemma [9] the function S determines 

- aky 



\ua{tj,Xk;,yj,r]j)\g\^/^{xk^)\ 



As Oj, i G Z+, are disjoint, this determines the index kj and thus also the 
point Xkj- Moreover 

hm Xfc, = hm -f{tj;yj,r]j) = 7(r(x, 0; a;, 0- 

□ 

The following result follows from Remark 2. 
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Lemma 10. Let us denote by X either SQ or d^SM. Let {yo,Vo) ^ 
satisfy 

riyo,Vo) < oo, iiT{yo,r]o);yo,r]o) ^ T^dM, 

where z = 7(r(?/o, ''7o); Vo, Vo)- Then there is a neighborhood V C X of {yo, tjq) 
such that T = i in V, where £ : U ^ I is the map of boundary intersection 
times defined in Remark^ for neighborhoods U <Z X and I C of {yo,T]o) 
and T{yQ,rio). In particular, r is smooth in V. 

Lemma 11. The set of (x, ^) such that 7(-; x,^) is transverse to DM is open 
and dense in 

dSM := e SM] X e DM}. 

Proof. As d-SM U d+SM is open and dense in dSM, it is enough to show 
that the set of (x,^) such that 7(-; x, ^) is transverse to dM is open and dense 
in d±SM. By the parametric transversahty theorem, see [20l Thm. 3.2.7], 
the claim follows from the fact that the evaluation map 

F^"" : d±SM xR^W 
F^^ ■.{x,U)^l{t-x,0 

is transverse to dM. □ 

Lemma 12. Let (a;o,^o) ^ -D(S). Then there is a sequence {{xj,^j))'^-^ C 
such that 'y{-;Xj,C,j) is transverse to dM and 

lim(xj,^j) = (xo,^o), hm T{xj,Q = r(xo,^o)- 

Proof. We denote tq := t{xo,C,o) and 

(zoXo) ■■= (7(^o;a;o,^o),-7(T"o;a;o,^o))- 

Remark [2] gives a map of boundary intersection times i : U I for neigh- 
borhoods U C SM."" and / C M of (2:0,^0) and tq. By Lemma 11 there is a 
sequence {{zj, Cj))j^i C SM fl U converging to {zq, Co) such that 7(-; zj, Q) is 
transverse to dM. 

We define tj := i{zjXj) and 

{xj,Q := {-f{tj;ZjXj),-i{tfZjXj))- 
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Then {xj,$,j) — )■ (a;o,^o) as j — )■ oo. In particular, there is J > 1 such that 
{xj,^j) e d-SM for j > J. By Lemma 

r(xo,^o) < hminf r(xj, ^j) < hmsup r(xj, ^j) 



< hm iizjXj) = ((zoXo) = t{xo,^o). 



□ 



Lemma 13. Let (xo,^o) G -D(S) be such that 7(-;a;o,^o) 'is transverse to 
dM. Then there is {yo,rio) G 5*^2 lying on the geodesic 7(-;xo,^o) o^nd a 
neighborhood V C Sy^Q of tiq such that the following conditions hold. 

(CI) The map rj i-> T{yQ,ri) is smooth V — (0,oo). 

(C2) The map 

{x{r]),^{r])) := {-f{T{yo,r]);yo,r]),'y{T{yo,r]);yo,r])) (22) 

is smooth V -D(S) and (a;(?7o), ^(^7o)) = (a^o,^o)- 
(C3) The map 

eiv)-=rixiv),av)) + riyo,v) (23) 
is smooth V — )■ (0, oo). 

(C4) There is a neighborhood W C dM o/7(r(xo, ^o); a^o, ^o) such that 

V^l{r{x{7]),^{v))]x{r]),^{r])) (24) 

is a difjeomorphism V ^ W. 



Proof. We denote 7(t) := 7(t;xo,^o) and Zq := 7(t(xo,^o))- By remark [2 
7(— t) G Q for small t > 0. Moreover, the points that are conjugate to zq 
along 7 are discrete on 7, see e.g. |24j . 
Thus there is tq > such that 

{yo,Vo) ■= (7(-^o),7(-^o)) 
is in 5*^2, yo is not conjugate to Zq along 7, t(?/o,''7o) = and 

(7(t"o; 1/0,^0), 7(^0; 1/0,^0)) = (a;o,^o)- 
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By Lemma 10 there is a neighborhood Vq C Sy^Vl of r]Q such that r] k-)- 
T{yo,ri) is smooth in Vq. Hence the function r] h-> (a;(?7), ^(r^)) maps r/o to 
{xo,^o) and is smooth in Vq. Moreover, this smoothess, tranversahty of 



li'^^o^^o) and Lemma 10 imply that there is a neighborhood Vi C Vq of 770 
such that {x{ri),^{ri)) G d^SM and i] 1— )■ r(x(?7), ^(r/)) is smooth Vi — > (0, 00). 
In particular, (^(r/), ^(77)) G -D(S) for all 77 G Vi. We have shown that (?/o, ^70) 
and 1/1 satisfy (C1)-(C3). 
We have 

{l{s;yo,v),i{s;yo,v))\s=t+r{yo,v) = {l{t;x{r]),C{r])),'y{t;x{7]),C{r]))). (25) 
In particular, 7(£(?7o); ?/o, ''7o) = and 

Moreover, as ?/o is not conjugate to zq along 7, there are neighborhoods 
^2 C Vi, Jo C (0,00) and Uq C M" of rjo, l{r]o) and 2:0 such that {t,r]) 
7(t; 7/0, v) is a diffeomorphism V2 x Jq — )■ Uq. 

There is a neighborhood C V2 of ?7o such that ^(V^) C /q. The graph of 
7] I— )■ £(77) is an (77 — 1) dimensional submanifold on x Jg. Hence the diffeo- 
morphism (t, rj) H- 7()f:; t/qj ''?) maps it onto a (77 — 1) dimensional submanifold 
W oHJq. Moreover, zq E W and W C dM. Thus 1^ is a neighborhood of zq 
in (9M. □ 

Lemma 14. Let (iCoj^o) ^ -^(S) ^^^'^ {yojVo) ^ ^Ti satisfy conditions (Cl)- 
(C4) of Lemma [igj /or neighborhoods V C a?T-c? ly C c?M 0/7/0 o^c? 



Zo := 7(T(a;o, ^0); 2:0, ^o)- We denote by F : W ^ V the inverse map of (24-) ■ 
Then 

grads,SoF)U^,^^=jl, (26) 



where £ : V ^ {0, 00) is the function (23) and 7^ is the orthogonal projection 
o/7(r(a;o,^o);xo,^o) intoT,^,dM. 

Proof. Let a : (— e, e) ^ W he a. smooth curve such that cr(0) = Zq. We 
define 

r : (-e, e) X M ^ T{s, t) := ^{t; 7/0, F{a{s))). 



We denote A := £ o F o a and io := £(?7o)- By equation (25) 

r(s,A(s)) = 7(r(a;(77),^(77));x(77),^(r/))|^=i.(^(,)) = (t(s), 
{dtT){0,io) = 7(^;2/o,^?o) = 7(7-(xo,^o);a;o,^o)• 
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Hence 

dr(O) = 9,r(s, A(s))|,=o = {dsTmio) + (W(0,4)A'(0). 
The curve t T{s, t) is a unit speed geodesic for all s G (— e, e). Hence 
(a(0),7(r(xo,eo);xo,eo))5 = ((5.r,W, + A'(0)(9,r,W,)|,^o^^^^^^ (27) 

We define 

C{s,l):= [ \dtT{s,t)\gdt, (s,/) G (-e,e) X (0,oo). 

Then £(s, 0; -s G (— e, e) is the length of a unit speed geodesic on the interval 
[0, /]. Thus jC(s, /) = / for all s G (— e,e). We may derive an expression for 
dsC{s,l)\s=o as in [36l Prop. 6.5] 



dsCis,l)l=o = I {DtdsT,dtT)gdtl=o. 
Jo 



As t H- r(s,t) is a geodesic, DtdtT{s,t) = and thus 

dt{d,r,dtr)g = {DtdsT,dtT)g. 

Moreover, r(s,0) = yo for all s G (— e,e) and thus 9sr(s,0) = 0. Hence 

= d,C{s,l)l=o= [ dt{dsT,dtT)gdt\s=o = {d,T,dtT)g\s=o,t=i, /g(0,cx)). 
Jo 

By ( [27) ) we have 

((^(0), tJo)^ = (fT(0), 7(r(xo, ^0); xo, ^o))g 

= A'(0) = (d(£oF)|,=,„,a(0)) 

for all smooth curves a inW such that cr(0) = Zq, which proves the claim. □ 
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Theorem 4. The functions r : d^SM — )■ (0, oo] and 

z:D{^)^dM, ;^(x,0:=7(r(x,0;a:,e) 
together with the Riemannian manifold {VL,g\Q) determine 

7(r(x,0;a;,e), eZ)(S). 

Proof. The functions r and z on -D(E) determine the set B of points (xq, ^o) ^ 
such that the conditions (C1)-(C4) of Lemma 13 hold for some (?/o, Vo) ^ 

SVt. 

Let (xo,Co) ^ -B. We denote Co := i{'^{'^o-,io)]Xoiio)- The map 

determines its local inverse. Hence r and z determine the function F of 
Lemma 
vector 



14 



and thus they determine 7^ by the formula (26). As Co is a unit 



Co = il + (1 



17; 



T |2\l/2 



where z/^p the unit exterior normal vector of dM. Hence r and z determine 



Co for aU (xc^o) e B. 

Let {xq,^q) e -D(S). By Lemmata 
B such that 



12 



and 



13 



there is a sequence {{xj,^j))'jLi C 



hm (xj, ^j) = {xo, ^0), lim r(xj, ^j) = t{xo, ^q). 

Moreover, the functions r and z determine the set of such sequences and thus 
they determine 

lim ^{T{xj,Q);Xj,Q) = ^{t{xo,^o);xo,^o)- 

□ 

Theorems [2| [3] and |4] prove Theorem [T] formulated in the introduction. 
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